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Abstract 

A quasi-particle theory for monatomic gases in equilibrium is formulated and evalu- 
ated to yield the exact virial contributions to the thermodynamic state functions in 
lowest order of the density. Van der Waals blocking has necessarily to be accounted 
for in occupation number statistics. The quasi-particle distribution function differs 
from the Wigner function by a bilinear functional thereof. The progress made so far 
is promising with respect to a corresponding version of kinetic theory. 
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1 Introduction and summary 



The quasi-particle concept serves as a camouflage of interaction in many-body systems. 
The aim is to explain a macroscopic system in equilibrium as an ensemble of noninteract- 
ing and countable units called " quasi-particles" . Non-equilibrium may then be described 
by a kinetic theory for these quasi-particles. 

In this article the interaction-free theory is occupation number statistics which yields the 
entropy as a functional of the quasi-particle distribution function /. Interaction between 
gas atoms is taken care of by modelling the elementary cell volume and by suitable 
constraints for /. 

Our theory reproduces the exact quantum mechanical density corrections to the distri- 
bution function and to equilibrium thermodynamics. The only partial success of related 
attempts (e.g. [1]) stems from an insufficient handling of the combinatorial entropy, i.e. 
the neglect of van der Waals blocking. 
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2 One-particle distribution function 



The one-particle distribution function f(p ± , f, t) of gas kinetics enables the calculation of 
densities (in position space) for macroscopic quantities by taking moments over momen- 
tum space. We insist the two simplest ones be the number density n 

n(r,t)= J d 3 Pl /(p^t) (1) 
and the density of linear momentum 7? 

7f(r,t) = /d 3 Pl Pi/(Pi,r,*) . (2) 

Any justification of gas kinetics within the frame of a more fundamental theory starts 
with a definition of / which one is free to choose, provided the interpretation of the 
moments (1) and (2) is correct. Our quantum statistical ansatz reads 

/(P!,r,t) := f w {Pi,r,t) + V{p 1 ,r,t) , (3) 
where fw is the Wigner distribution function 

fwfa, r, t) = (2nhy 3 J dV ( ^ f (r + 1 t)^(r - T -, t) ) exp {l p x ■ rj (4) 
and \1/ is a bilinear functional thereof 

*(Pi,r,t):= (5) 

J d 3 p 2 dVidV 2 E (P, P)$ (P ~ P) [fw(Pi,r, t)f w (p 2 , r, i) - f w (p[, r, t)f w (p 2 , r, t) 

+ J d 3 p 2 dVidV 2 0(p, p)5 (p - P ) [f w (pi, r, t)fw(p 2 , r, t) + fw(p[,r, t)f w (p 2 , r, f)] , 

with P, P and p, p denoting centre-of-mass and relative momenta throughout the article 

P = Pi+P2 , P= \ (P1-P2) • (6) 
The even kernel 

E(p, p) = E(p, p) = E(-p, -p) (7) 
and the odd kernel 
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0(p,p) = -0(p,p) = 0(-p, -p) 



(8) 



remain to be specified. The structure of \I/ reminds of Boltzmann's collision integral 
although there is no 5-function for the kinetic energies of relative motion E p , E p > which 
are defined by 



->2 ->2 

E P = ^ = P ~. (9) 

2m re i m 

In any case equations (1) and (2) hold: 

\t) = (^(r,t)iP(r,t)) = J d 3 Pl f(p 1: r,t) = J d 3 Pl f w i$^t) (10) 



n[r, 



and 



n(r,t) = J d 3 Pl Pi /(Pi,r,t) = J d 3 Pl p x f w i$ x ,r,t) . (11) 

There is not only a quantum statistical motivation to the definition (3), but it also allows 
a quasi-particle interpretation. 



3 Quantum statistical background 



Explicit expressions for the kernels E and O follow from the theory of Kadanoff and 
Baym [2] , from which in the spatially homogeneous case and in T-matrix approximation 
the kinetic equation 



dt{fw + ^) = J B [fw\ (12) 

is obtained [3,4] with Jsifw] as Boltzmann's collision integral depending on f w . Here \I/ 
is just the functional (5) with the kernels 



0(p,p') = ^(2nh) 3 V(E p -E p ,) (13) 



{\(p\T ± (E p , + ie) \p>)f-\(p>\T ± (E p + i e) \p)\ 2 } 



and 
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E(p,p') = n(2nhf5(E p -E pl ) (14) 
Im((p| T ± (E p , + ie) \p')*(p\ T^E^+ie) \p')) 



+ ^(2irh) 3 V'(E p - E p , 



{\(p\T ± (E p , + ie) \p>)f + \(p>\T ± (E p + ze) |p)| 2 } . 

The T-matrix occurring here is the properly symmetrized momentum representation of 
the two-particle operator 



T(z) = V- Vj^—V , r\z) = ±- T(z) , (15) 

with H = Hkm + V being the Hamiltonian of relative motion. V is the principle value 
distribution and V it's derivative. The upper sign in T± (and elsewhere) refers to bosons 
and the lower sign to fermions. 

In equilibrium, quantum statistical mechanics yields the density expansion for / = f w +ty, 
which up to second order reads 



UA) = (27rm " r) 3/2 e " A (! + < 2B ^) + ^l))) ' ( 16 ) 

with k being Boltzmann's constant and 



4>(p 1 ) = ±\\T) e-^kr - A 3 (T) J d 3 p 2 e -2ir|^ + G(p)} . 
Here we have introduced the thermal wave length 



A = ^^. (17) 

The quantities F and G as well as the second virial coefficient B(T) are given as func- 
tional of the T-matrix: 



F(p)= l -{2KhfRe{(p\T ± {E p + ie) \p)) , (18) 



G( P ) 



(2nh) 3 J d 3 q 8{E p -E q ) 

*\m({p\T ± {E q + ie) \q)(q\ T±(E q + ie) \p)*) 



(19) 
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B(T) = B (T) + B 1 (T) + B 2 (T) (20) 
B (T) = T 2- 5/2 X 3 , B 1 (T)= { -^p-, B 2 (T) = ((G)), 
where (( • )) denotes the thermal average, e.g. 



((f)) = j p 3 P1 . (21) 

Jd v exp{-^} 

These formulae are exact if the two-particle interaction does not allow any bound states. 
Bi essentially accounts for long-range attraction and B 2 for hard repulsion, this corre- 
spondence being most striking in the van der Waals limit (cf. section 6). 



4 Quasi-particle picture for equilibrium 

We just work out the usual idea: 

(1) The entropy density s is represented as a functional of the one-particle distribution 
function. This is an outcome of occupation number statistics (combinatorial entropy) 

s = -kJ d 3 Pl [ /&) In (^(Pi)/(Pi)) (22) 

=F ± In (l ± ^(Pi)/(Pi)) 

with v e i as the volume of an elementary cell in six-dimensional ^u-space, v e i accom- 
modating one single-particle quantum state. Eq. (22) is well known as a standard 
result for non-interacting particles. However, the choice of v e \ and the constraints for 
/ may provide a camouflage of interaction. 

(2) For equilibrium the distribution function / is the one which minimizes s subject to 
appropriate constraints. 

(3) The constraints are a given number density n, 

n = J d 3 Pl ffa) , (23) 

and a given energy density u, 

u = J d 3 Pl efo) ffa) . (24) 

The quasi-particle interpretation is now introduced by way of ansatz (eqs. (25), (26)), its 
aim being to account for interaction effects in lowest order of the density. Strong repulsion 
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reduces the freely accessible volume for gas particles. Because of this effect, called "van 
der Waals blocking", more than just (2iTh) 3 is needed as an elementary cell volume 



(p 1 ) = (2 7 rn) 3 [l + Jd 3 p 2 G(\p\) f(p 2 ) 



(25) 



Also, an interacting gas particle carries with it a correlation cloud giving rise to an 
interactive contribution which changes the kinetic energy of a particle into the energy of 
a quasi particle 



For our variational problem the functions G and F are considered as given though, for the 
time being, unknown. They will be determined afterwards by comparing the equilibrium 
solution / = / e g with the corresponding expression from many-body quantum theory. 
The resulting thermodynamics then serves as a further touchstone of the quasi-particle 
interpretation. 

The solution of the variational problem is obviously equivalent to 



with $ and a as Lagrange parameters due to the the constraints for /. Comparison with 
the thermodynamic identity for the entropy S at constant volume 



reveals that $ means the temperature, ■& — T, and a means the chemical potential, a = //. 
According to (27) / = f eq is equivalently determined by the fixed-point equation 




(26) 




(27) 



dS = I (dU - /idJV) 



(28) 




(29) 



with 




C = exp {fin} 



(30) 



and 
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K{p 1 ) = J d 3 p 2 ^2f(p 2 )F(\p\) (31) 



5 Lowest-order density corrections 

The fugacity ( can be given as a power series in n 

( = n\ 3 (l + 2nB(T) + ---) (32) 

and the fixed point equation (29) may be iterated starting off with a Maxwellian nor- 
malized to n. This yields a density expansion according to which - apart from third and 
higher order contributions - one regains eq. (16) for /, but now with 

= ±A 3 e" ^ (33) 

\3 g 

"(2^F Jd 3 p 2 e-^{PF(\p\) + G(\p\)} ■ 
Therefore, the quasi-particle picture independently explains eq. (16), if and only if 

F(p) = F(p) (34) 

and 

G{p) = G(p) (35) 

Having determined / up to second order in n, we deduce the entropy with its lowest order 
density corrections from eq. (22): 

s = riK 1 1 - In (nA 3 ) - n [B(T) + TB'(T)} J , (36) 

with B(T) defined by eqs. (20). Analogously evaluating eq. (24), we obtain the two leading 
contributions to the energy density: 

u = ukT jJ-nT^T)} . (37) 
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The last two results imply the pressure equation of state 



( du \ 

p = ni— J - u = ukT {I + nB(T)} . (38) 

\ / s/n 

The exact density corrections, i.e. virial contributions have thus been obtained. 



6 Classical van der Waals approximation 



Considering distinguishable particles, one has to neglect quantum statistical contribu- 
tions. This means 



B = . (39) 
Then with (20), (34) and (35) 



B(T) = y^p- + {{G)) (40) 

holds and suggests a comparison with the van der Waals version of the second virial 
coefficient 



B vdW (T) = --^ + b, (41) 
which is readily obtained from the model equation of state 

(p + n 2 aj (1 - rib) = ukT (42) 
and the corresponding density expansion 



p = ukT ( 1 + n 



kT. 



+ o(n 3 ) . (43) 



The van der Waals limit therefore obviously means 



(( F )) = const = -o and ((G)) = const = b . (44) 

This model assumption is actually quite reasonable as we are going to demonstrate for 
4 He atoms interacting via a Lennard- Jones-potential lacking bound states [5]: 
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V(r)=4V 



(T\ 12 /CT\ 

r / 



V = 10.22 K k 



2.56 A . 



(45) 



Because the interaction is radially symmetric the following relations between the (anti-) 
symmetrized T-matrix, the scattering amplitude f±(p,9) and phase shifts 5i(p) 



f±{P-,Q) = -K 2 mh{p\ T±(E p + ie) \q) , \p\ = |g| , p ■ q = \p\\q\cos(6) 

/±(p, *) = -£' (2/ + 1) e t5lip) sin(<S,(p)) P(cos(0)) 
Pi 



(46) 
(47) 



may be used (see e.g. [6]), where the summation runs over even I for bosons and odd / for 
fermions. F(p) and G(p) (compare to eqs. (18) and (19)) can then be expressed in terms 
of phase shifts which is compatible with the Beth Uhlenbeck result for B(T) [7] 



F(p) 
G(P) 
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Fig. 1. F and G as functions of ka (thick solid lines) and thermal weight functions 
w oc p 2 exp {— Ep/ (kT)} for two temperatures. 

Figure 1 nicely shows that for a large region of temperatures ((F)) and ((G)) may be 
nearly considered as constants, i.e. not depending on temperature. In this case the single 
particle energy, the elementary cell volume and the mean energy are easily determined as 



->2 



{2nhf[l + bn} , u = nnT^-^ 



(49) 
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7 Physical meaning of the quasi-particles 

One may imagine a quasi-particle to be a gas particle together with its surrounding 
correlation cloud which is described by the radial distribution function g(r,T). This in- 
terpretation suggests itself because the virial correction to the energy of the ideal gas (eq. 
(37)) is mainly determined by g(r, T) in an obvious way. The interpretation is immediately 
evident in the classical case where 




(50) 



with 




(51) 



as the first term of the density expansion 



g(r,T)=g (r,T) +ng 1 (r,T) +n 2 g 2 (r,T) H . 



(52) 



Therefore 




(53) 



is the classical virial correction to the internal energy. 



In general, however, for a homogeneous system 



g(r, T) = \{{ ^{r') ^(r' + r) ^(r + r) V>(r) )) 



(54) 




(55) 



with ff ree being the free energy per particle, whence 



g (r,T) = 2KT 



5B(T) 
SV(r) 



(56) 



Then quantum mechanically [8,4] 




(57) 
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where the scattering eigenstates \p)( ) satisfy the Lippmann Schwinger equation which 
reads in momentum representation 

( $ | p ) H = 8(p' - p) - (E p , -E p - ie) ( p' | V | ?)<->. (58) 
The virial correction to the energy density (eqs. (24), (37))is made up of three terms 

n\T TB'iT) = J d 3 Pl |^*m(?i) + / d'p, d 3 p 2 / M ( Pl ) A,(p 2 ) F(p) (59) 

+ /dV^fe(?,). 

Here /m is the Maxwellian normalized to n and ^ m is our functional (5) with f w replaced 
by fu- fw,2 denotes the second term in the density expansion of the Wigner function 

fw(Pi) = Jm(Pi) + n2 fw,2(Pi) + ■■■ ■ (60) 
Looking more closely at eq. (59) and taking account of (57) one can show that 

I d3pi 2^ MiS>l) + / d " Pl d " P2 MA) /m&) = / T)V{r) ■ (61) 

Therefore, in the classical limit the last term in eq. (59) must vanish. 
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